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Exploiting quantum properties to outperform 
classical ways of information-processing is an out¬ 
standing goal of modern physics. A promising 
route is quantum simulation, which aims at imple¬ 
menting relevant and computationally hard prob¬ 
lems in controllable quantum systems. Here we 
demonstrate that in a trapped ion setup, with 
present day technology, it is possible to realize a 
spin model of the Mattis type that exhibits spin 
glass phases. Remarkably, our method produces 
the glassy behavior without the need for any dis¬ 
order potential, just by controlling the detuning 
of the spin-phonon coupling. Applying a trans¬ 
verse field, the system can be used to benchmark 
quantum annealing strategies which aim at reach¬ 
ing the ground state of the spin glass starting 
from the paramagnetic phase. In the vicinity of a 
phonon resonance, the problem maps onto num¬ 
ber partitioning, and instances which are difficult 
to address classically can be implemented. 


Spin models are paradigms of multidisciplinary science: 
They find several applications in various fields of physics, 
from condensed matter to high energy physics, but also 
beyond the physical sciences. In neuroscience, the fa¬ 
mous Hopfield model describes brain functions such as 
associative memory by an interacting spin system pj. 
This directly relates to computer and information sci¬ 
ences, where pattern recognition or error-free coding can 
be achieved using spin models Q. Importantly, many 
optimization problems, like number partitioning or the 
travelling salesman problem, belonging to the class of 
NP-hard problems, can be mapped onto the problem of 
finding the ground state of a specific spin model SI- 
This implies that solving spin models is a task for which 
no general efficient classical algorithm is known to ex¬ 
ist. In physics, analytic replica methods have been de¬ 
veloped in the context of spin glasses @, 0. A contro¬ 
versial development, supposed to provide also an exact 
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numerical understanding of spin glasses, regards the D- 
Wave machine. Recently introduced on the market, this 
device solves classical spin glass models, but the under¬ 
lying mechanisms are not clear, and it remains an open 
question whether it provides a speed-up over the best 
classical algorithms 0-0]. 

This triggers interest in alternative quantum systems 
designed to solve general spin models via quantum sim¬ 
ulation. A noteworthy physical system for this goal are 
trapped ions: Nowadays, spin systems of trapped ions 
are available in many laboratories fiol-llfij , and adiabatic 
state preparation, similar to quantum annealing, is ex¬ 
perimental state-of-art. Moreover, such system can ex¬ 
hibit long-range spin-spin interactions [161 ] mediated by 
phonon modes, leading to a highly connected spin model. 
Here we demonstrate how to profit from these properties, 
using trapped ions as a quantum annealer of a classical 
spin glass model. 

We consider a setup described by a time-dependent 
Dicke-like model: 


Ho (t) = ^2 Tv^malnClm + ^ Mlj 
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with a m annihilating a phonon in mode m with frequency 
ujm and characterized by the normalized collective coor¬ 
dinates . The second term in Hq couples the motion of 
the ions to an internal degree of freedom (spin) through a 
Raman beam which induces a spin flip on site i, described 
by <t* , and (de)excites a phonon in mode m. Here, fl,; 
is the Rabi frequency, &u rec oil the recoil energy, and wl 
the beatnote frequency of the Raman lasers. Before also 
studying the full model, we consider a much simpler ef¬ 
fective Hamiltonian, derived from Eq. (ID by integrating 
out the phonons jl6l - Il9| . The model then reduces to a 
time-independent Ising-type spin Hamiltonian 


H J = -hY J J U°i4- ( 2 ) 

ij 


Each phonon mode contributes to the effective coupling 
Jij in a factorizable way, proportional to £im£jm, and 
weighted by the inverse of the detuning from the mode 
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The £i m imprint a pattern to the spin configuration, sim¬ 
ilar to the associative memory in a neural network 111I20l|. 
The connection between multi-mode Dicke models with 
random couplings, the Hopfield model, and spin glas s 
physics has been the subject of recent research [21 (231]. 
and the possibility of addressing number partitioning was 
mentioned in this context [23 ]. 

Before proceeding, we remind the reader that the con¬ 
cept of a spin glass used in the literature may have 
different and controversial meanings (cf. [2J-|27(): i) 
long-range spin glass models [b] and neural networks 
28, [29], believed to be captured by the Parisi picture 
30, [Ml and replica symmetry breaking. These lead to 
hierarchical organization of the exponentially degener¬ 
ated free energy landscape, breakdown of ergodicity and 
aging, slow dynamics due to a continuous splitting of 
the metastable states with decreasing temperatures (cf. 
fs3]). ii) Short-range spin glass models believed to be 
captured by the Fisher-Huse [33] droplet/scaling model 
with a single ground state (up to a total spin flip), but a 
complex structure of the domain walls. For these models, 
aging, rejuvenation and memory, if any, have different 
nature and occurrence [H, [H[; iii) Mattis glasses [35} . 
where the huge ground state degeneracy becomes an ex¬ 
ponential quasi-degeneracy, for which finding the ground 
state becomes computationally hard ([36|, see Subsec¬ 
tion “Increasing complexity”). Note that exponential 
(quasi)degeneracy of the ground states (or the free en¬ 
ergy minima) characterizes also other interesting states: 
certain kinds of spin liquids or spin ice, etc. 

Here we analyse the trapped ion setup. Even with¬ 
out explicit randomness (Hj = Q = const.), the coupling 
to a large number of phonon modes suggests the pres¬ 
ence of glassy behavior. This intuition comes from the 
fact that the associative memory of the related Hopfield 
model works correctly if the number of patterns is at 
most 0.138IV, with N the number of spins (28|. For a 
larger number of patterns, the Hopfield model exhibits 
glassy behavior since many patterns have similar energy 
and the dynamics gets stuck in local minima. However, 
it is not clear a priori how the weighting of each pat¬ 
tern, present in Eq. ©, modifies the behavior of the spin 
model. In certain regimes the detuning suggests to ne¬ 
glect the contributions from all but one mode, leading 
to a Mattis-type model with factorizable couplings [35[ . 
Jij oc £im£jm- Strikingly, the possibility of negative 
detuning, i.e., antiferromagnetic coupling to a pattern, 
drives such system into a glassy phase, characterized by 
a huge low-energy Hilbert space. The antiferromagnetic 
Mattis model is directly connected to the potentially NP- 
hard task of number partitioning j3§,[37j- Its solution can 
then be found via quantum annealing, i.e. via adiabatic 
ramping down of a transverse magnetic field, see also the 
proposal of Ref. [38| for a frustrated ion chain. 


(a) S = -2tt x 159 kHz (5) S = 2wx 796 kHz 



FIG. 1: Power-law decay vs. quasi-random patterns. 

We plot the coupling constants Jij (in units of J rma = 
n(n-i) l J o| 2 )> I° r a system of six ions at different de¬ 

tunings (5 = u>n -wl. Rabi frequencies are taken as constants 
Qi = const, (a) For negative detuning, 8 — 2i r x 159 kHz, 
interactions have a power-law decay, (b) For positive detun¬ 
ing, 8 = 2n x 796 kHz, the coupling constants resemble a spin 
glass. 


We start by giving analytical arguments to demon¬ 
strate glassy behavior in the classical spin chains. Us¬ 
ing exact numerics, we then focus on the quantum Mat¬ 
tis model. By calculating the magnetic susceptibility, 
and an Edward-Anderson-like order parameter, we dis¬ 
tinguish between glassy, paramagnetic, and ferromag¬ 
netic regimes. The annealing dynamics is investigated 
using exact numerics and a semi-classical approximation. 
We demonstrate the feasibility of annealing for glassy in¬ 
stances. Finally, we show that the memory in the quan¬ 
tum Hopfield model is real-valued rather than binary. 
This property might be useful for various applications 
of quantum neural networks such as pattern recognition 
schemes. 


Results 


Tunable spin-spin interactions 


We start our analysis by inspecting the phonon modes. 
For a sufficiently anisotropic trapping potential, the ions 
self-arrange in a line along the z-axis j39j, (4Q]. The 
phonon Hamiltonian H p h is obtained by a second-order 
expansion of Coulomb and trap potentials around these 
equilibrium positions zp. H p h = {m/2) with 

qi the displacement of the ith ion from equilibrium. For 
the transverse phonon branch, is given by fl6j 
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Our exact numerical simulations have been performed for 
six 40 Ca + ions in a trap of frequency co± = 2n x 2.655 
MHz, as used in a recent experiment (l4j . In order to 
maximize the bandwidth of the phonon spectrum and 
thereby facilitate the annealing process, we choose the 
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radial trap frequency uj z as large as allowed to avoid zig¬ 
zag transitions, uj z < 1.37a; j . • N~ 0 - 86 . Diagonalizing the 
matrix Vij leads to the previously introduced mode vec¬ 
tors = (£1 m, ■ ■ ■, £,Nm ), which are normalized to one, 

and ordered according to their frequency a; m : 

= (5) 

The mode with largest frequency, a in = a>j_, is the 
center-of-mass mode, = N _1/,2 . Parity symmetry of 
Vij is reflected by the modes, £,;- m = ±£(jv+i_i) m , an d 
even (+) and odd (—) modes alternate in the phonon 
spectrum. We focus on even N, for which all components 
£ i m are non-zero. Except for the center-of-mass mode, all 
modes fulfill = 0. 

Previous experiments flQj - lp| have mostly been per¬ 
formed with a beatnote frequency u l several kHz above 
wjv, leading to an antiferromagnetic coupling Jij < 0 
with power-law decay, see Fig. QJa). Despite the pres¬ 
ence of many modes, the couplings J,; 7 then take an or¬ 
dered structure. This work, in contrast, focuses on the 
regime wl < wjv, where modes with both positive and 
negative Sj lrn contribute, and ferro- and antiferromag¬ 
netic couplings coexist, cf. Fig. [jjb). This reminds of 
the disordered scenario of common spin glass models like 
the Sherrington-Kirkpatrick model [5j. In the following 
we study the properties of the time-independent effec¬ 
tive spin model, before considering the full time-depedent 
problem involving spin-phonon coupling. 

Classical Mattis model 

Close to a resonance with one phonon mode, simple 
arguments allow for deducing the spin configurations of 
the ground states. In this limit, we can neglect the other 
modes, leading to a Mattis model. A single pattern £ m 
then determines the coupling, J^ oc £,im^,jm- The sign of 
Jij depends on the sign of the detuning: Below the reso¬ 
nance, we have sign(Jy) = sign(^ m ^j m ), and accordingly 
the energy —hJija x a J x is minimized if a l x and a 3 . are ei¬ 
ther both aligned or both anti-aligned with £,; m and £j m . 
Thus, we have a two-fold degenerate ground state given 
by the patterns ±[sign(£i m ),..., sign(£jy m )]. We refer to 
this scenario as the ferromagnetic side of a resonance. 

Crossing the resonance, the overall sign of the Hamil¬ 
tonian changes. Naively, one might assume that this 
should not qualitatively affect the physics, since we have 

= -jf -> o, that is, there is an equal 
balance between positive and negative Jij. This expec¬ 
tation, however, turns out to be false. Recalling the re¬ 
lation between the Mattis model and number partition¬ 
ing Sana, the antiferromagnetic model maps onto 
an optimization problem in which the task is to find the 
optimal bi-partition of a given sequence of numbers (£,),, 

such that the cost function E = £* — X^'et £j) 

minimized. Here, the two partitions are denoted by f and 
4-. For a Hamiltonian of the form H = J^ij ^i a x^j a i = 


(Y^i^i a x ) 2 1 eigenvectors of a x are Hamiltonian eigen¬ 
states with an energy precisely given by the cost function 
E. Thus, in the limit of just one antiferromagnetic reso¬ 
nance, the ground state of H is exactly the configuration 
that minimizes the cost function. 

With this insight, the ground states of the spin model 
are easily found exploiting the system’s parity symmetry: 
For even modes, = £(jv+i_i) m , and we simply have to 
choose (<r®) = — (cr x +1 ~' 1 ) to minimize the cost function. 
For odd modes, = —£(jv+i-i)mi and we must choose 
Wx) = ( &x +1 ~*)■ I n both cases, this implies that we can 
choose half of the spins arbitrarily, leading to at least 
2 n / 2 ground states. 

The important observation is that an exponentially 
large number of ground states exists in the limit of being 
arbitrarily close above a resonance. This is a characteris¬ 
tic feature of spin glasses, yet it does not lead to computa¬ 
tional hardness. In fact, as pointed out in Ref. [221], the 
number partitioning problem with exponentially many 
perfect partitions belongs to the “easy phase”. How to 
reach hard instances will be explained in the section be¬ 
low. 

Pushing our arguments further we consider the influ¬ 
ence of a second resonance: In between two resonances, 
the exponential degeneracy of the antiferromagnetic cou¬ 
pling on one side is lifted by the influence of the ferromag- 
netically coupled mode on the other side. Interestingly, 
this does not lead to frustration, since even- and odd- 
parity modes alternate in the phonon spectrum, and the 
pattern favored by the ferromagnetic coupling is always 
contained in the ground state manifold of the antiferro¬ 
magnetic coupling. Accordingly, between two modes the 
ground state pattern is uniquely defined by the upper 
mode. 

Beyond this two-mode approximation, we rely on nu¬ 
merical results. Taking into account all modes, exact 
diagonalization of a small system (N < 10) shows that 
the two-mode model captures the behavior correctly: At 
any detuning, the degeneracy due to the nearest antifer¬ 
romagnetic coupling is lifted in favor of the pattern of 
the next ferromagnetic coupling. In Fig. [24a), we plot 
the cumulative density of states Pcum(-E'), that is, the 
number of states with an energy below E. The corre¬ 
sponding phonon resonances uj m are marked in Fig. [24 b). 
The curves clearly reflect the very different behavior in 
the red- and blue-detuned regimes: Fig. [2] illustrates 
the quick increase of p C um{E!) at low energies, when the 
laser detuning is chosen on the antiferromagnetic side 
of a phonon resonance. In contrast, a low density of 
states characterizes the system on the ferromagnetic side 
of a resonance. In intermediate regimes, as shown for 
S = 2tt x 199 kHz, the spectrum is symmetric. 

A breakdown of the two-mode approximation is ex¬ 
pected for large numbers of spins: Since the distance be¬ 
tween neighboring resonances approximately scales with 
1/N (at fixed trap frequencies), the influence of addi¬ 
tional modes grows with the system size. The combined 
contribution of all antiferromagnetically coupled modes 
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FIG. 2: Cumulative density of states, (a) For six ions 
and at trap frequency co z = 2n x 770 kHz, we plot the num¬ 
ber of states (divided by the total Hilbert space dimension) 
below the normalized energy treshold (0: ground state en¬ 
ergy, 1: energy of highest state), for different detunings from 
the center-of-mass mode. The density of states at low ener¬ 
gies is seen to strongly increase when the detuning is slightly 
above a phonon resonance. In (b), the position of the reso¬ 
nances (measured from the center-of-mass mode at <5 m = 0) 
are shown. 


tries to select the fully polarized configurations as the 
true ground state, while all ferromagnetically coupled 
modes, except for the center-of-mass mode, favor fully 
unpolarized configurations. As a consequence, it is a pri¬ 
ori unclear which pattern will be selected in the presence 
of many modes. 

This observation is crucial from a point of view of com¬ 
plexity theory. In the presence of parity symmetry nei¬ 
ther the one-mode problem (i. e. the number partition¬ 
ing problem), nor the two-mode approximation are hard 
problems, as they can be solved by simple analytic argu¬ 
ments. However, when many modes lift the degeneracy of 
the exponentially large low-energy manifold in an a priori 
unknown way, one faces the situation where a true but 
unknown ground state is separated only by a very small 
gap. Identifying this state then usually requires scan¬ 
ning an exponentially large number of low-energy states, 
and classical annealing algorithms can easily get stuck in 
a wrong minimum. Below, we discuss how a transverse 
magnetic field opens up a way of finding the ground state 
via quantum annealing. Moreover, we will discuss strate¬ 
gies to make also the one-mode model, i.e. the number 
partitioning problem, computationally complex. 


Increasing complexity 

As discussed above, the instances of the number par¬ 
titioning problem realized in the ion chain are simple to 
solve due to parity symmetry. This is a convenient fea¬ 
ture when testing the correct functioning of the quantum 
simulation, but our goal is the implementation of compu¬ 


tationally complex and selectable instances of the prob¬ 
lem in the device. One strategy is the use of microtraps 
to hold the ions [4lj| . The equilibrium positions of the 
ions can then be chosen at will, opening up the possibil¬ 
ity to control the values of the £j m . The computational 
complexity of the number partitioning problem then de¬ 
pends on the precision with which the can be tuned. 
If the number of digits can be scaled with the number of 
spins, one enters the regime where number partitioning 
is proven to be NP-hard [371] - Thus, the number of digits 
must at least be of order log 10 N , which poses no problem 
for realistic systems involving tens of ions. 

Another way of enhancing complexity even within a 
parity-symmetric trap would be to “deactivate” some 
spins by a fast pump laser. For example, if all spins 
on the left half of the chain are forced to oscillate, 
crj —» cr^e l “ ;pumpt , the part of the Hamiltonian which 
remains time-independent poses a number partitioning 
problem of N /2 different numbers. 

Another promising approach was recently suggested in 
Ref. iHJ: Operating on the antiferromagnetic side of 
the center-of-mass resonance, single-site addressing al¬ 
lows one to use the Rabi frequency for defining the in¬ 
stance of the number-partitioning problem. This could 
indeed be the step to turn the trapped ions setup into 
a universal number-partitioning solver, where arbitrary 
user-defined instances can be implemented. 

If one is not interested in the number-partitioning 
problem itself, one might also increase the system’s com¬ 
plexity via resonant coupling to more than one mode. 
Equipping the Raman laser with several beatnote fre¬ 
quencies co and Rabi frequencies H, it is possible to 
engineer couplings of the form ItI : 


Mmax N 

Ja « E E 


‘V/mC, 


jm 


M=1 


LOm ~ U 
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( 6 ) 


With an appropriate choice of Rabi frequencies and de¬ 
tunings, this allows for realizing the Hopfield model, 
Jij OC ^im^jm.^1 where each coupling p is assumed 

to be in resonance with one mode m M . For ferromagnetic 
couplings, the low-energy states again are determined by 
the signs of the £j m , but in general the different low- 
energy patterns are not degenerate, and a glassy regime 
is expected for large ^ max [281 ]. 


Quantum phases 

So far, we have considered classical spin chains lacking 
any non-commuting terms in the Hamiltonian. Quantum 
properties come into play if we either add an additional 
coupling a y a y’ or a transverse magnetic field: 

H B = hBj2°l (?) 

i 

The latter has been realized in several experiments 1 1 (1 
12 [3 > an d is convenient for our purposes, as the field 
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B (Jrms) 


FIG. 3: Magnetic phases. Upon varying the laser de¬ 
tuning 5 and the transverse magnetic field B, we mark, for 
N = 6 ions and lo z = 2n x 770 kHz, those regions in configura¬ 
tion space where the ferromagnetic order parameter i'fm, the 
spin glass order parameter ^sg, or the longitudinal magnetic 
susceptibility x take larger than average values (Vfm = 0.13, 
zzsg = 5, x = 1-6). Below each phonon resonance (marked by 
the dashed horizontal lines), there is a regime where large val¬ 
ues of v sg and % indicate spin glass behavior for sufficiently 
weak transverse field B. In order to break the Z 2 symme¬ 
try of the Hamiltonian, all quantities were calculated in the 
presence of a biasing field ecrj. (with e = — J rm s). 


strength B , if decaying with time, provides an anneal¬ 
ing parameter: For large R, all spins are polarized along 
the ^-direction, whereas for vanishing B one obtains the 
ground state of the classical Ising chain. As argued 
above, the latter exhibits spin glass phases with an expo¬ 
nentially large low-energy subspace. Even in those cases 
where the true ground state is known theoretically, find¬ 
ing it experimentally remains a difficult task. Our system 
hence provides an ideal test ground for experimenting 
with different annealing strategies. 

Before presenting results for the simulated quantum 
annealing, let us first discuss the different phases ex¬ 
pected for the effective Hamiltonian H e g = Hj + Hb + 
ea ,The last term is a bias introduced to break the Zi 
symmetry. Our distinction between phases is based on 
certain quantities which combine thermal and quantum 
averages 


w*\ = Ea( a I • |A) a cxp(-A A /fc B r) 

; >T ~ Ex^M-Ex/k B T) ’ W 

with | A) denoting Hamiltonian eigenstates at energy Ex ■ 
For a = 1, ((')“)t reduces to the normal thermal av¬ 
erage. We will use low, but non-zero temperatures T 
of the order of the coupling constant, accounting in this 
way for the huge quasi-degeneracy in the glassy regime. 
The thermal average (-)t plays a role somewhat similar 
to the disorder average, as it averages over various quasi¬ 
ground states (pure thermodynamic phases). We there¬ 
fore expect ((<r* ))r to go to zero in the glassy phase. In 


contrast, a non-zero average ((<t*))t detects the ferro¬ 
magnetic phase of the Mattis model, while it vanishes in 
the paramagnetic state. Taking its square to get rid of 
the sign, we obtain a global ferromagnetic order param¬ 
eter by summing over all spins: 

J'fm = y^((c r s))r- (9) 

i 

On the other hand, in the spirit of an Edwards-Anderson- 
like parameter, we consider thermal averages of squared 
quantum averages, i.e. 


J'ea = YAK) 2 )t. (10) 

i 

At sufficiently low temperature this average would still 
be zero for a paramagnetic system, but now it remains 
non-zero for both ferromagnetic and glassy systems. Ac¬ 
cordingly, a parameter which is “large” only for glassy 
systems is given by the ratio i/sg = ^rS used in Fig. [3] 
to detect glassy regions. 

Thermal averages are difficult to measure, but the con¬ 
tained information is also present in linear response func¬ 
tions at zero temperature. Therefore, we have calculated 
the longitudinal magnetic susceptibility %, i.e. the re¬ 
sponse of the system to a small local magnetic field h J x 
along the o x direction (in units J rms )' 



Due to the (quasi)degeneracy in the glass, one expects a 
huge response even from a weak field, and thus a diver¬ 
gent susceptibility. 

We have calculated ffm, J'sGi and x, for N = 6 be¬ 
tween 0 < 2n x S < 8 MHz, and 0 < B < 4J rms . For 
the thermal averaging, we have chosen a temperature 
k B T = J rms . The results are summarized in Fig. [3l in¬ 
dicating the regions where these quantities take larger 
values than their configurational averages z/fm, I'sg, and 
X, defined as / = J dB f d Sf(B,S)/(B max S max ). In this 
way, we identify and distinguish ferromagnetic behavior 
below, and glassy behavior above each resonance. Re¬ 
gions of large susceptibility x overlap with regions of large 
v sg, attaining numerical values which are three orders of 
magnitude larger than the corresponding averages. For 
sufficiently strong field B , in contrast, none of these pa¬ 
rameters is large, indicating paramagnetic behavior. 

Note that the existence of the purported glassy phase 
in the quantum case is an open problem. We provide 
here the evidence only for small systems, since it is nu¬ 
merically feasible and corresponds directly to current or 
near-future experiments. If we increase the complexity of 
our system by resonant coupling to many phonon modes, 
as discussed in the last paragraph of the previous subsec¬ 
tion, the glassy behavior will result from the interplay of 
contributions of many modes similarly as in the Hopfield 
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model with Hebbian rule and random memory patterns. 
Here the beautiful results by Strack and Sachdev [21] 
the quantum analog of the Amit et al. i!29j machinery 
- can be applied directly to obtain the phase diagram 
for large N. If, however, we increase the complexity by 
random positioning of the ion traps, then the resonance 
condition will pick up the contribution from one domi¬ 
nant (random) mode, and the Hebbian picture will apply. 


Simulated quantum annealing 


We will now turn to the more realistic description of 
the system in terms of a time-dependent Dicke model, 
described by the Hamiltonian H$ in Eq. © with an addi¬ 
tional transverse field Hs{t) from Eq. ©. We assume an 
exponential annealing protocol B(t) = H max exp(— t/r). 
Again we apply a bias field /ibias = Secrj,, lifting the 
Z 2 degeneracy of the classical ground state. We study 
the exact time evolution under the Hamiltonian H(t) = 
Hp(t) + Hs(t) + hbia S , using a Krylov subspace method 
[4J], and truncating the phonon number to a maximum 
of two phonons per mode. 

Initially, the system is cooled to the motional ground 
state, and spins are polarized along er z . Choosing 
B max 3> e, Jrmsj this configuration is close to the ground 
state of the effective model Hj + Hb + /ibias at t = 0. If 
the decay of B(t) is slow enough, and if the entanglement 
between spins and phonons remains sufficiently low, the 
system stays close to the ground state for all times, and 
finally reaches the ground state of Hj. 

We have simulated this process for TV = 6 ions, as 
shown in Fig. H As a result of the annealing, we are 
not interested in the final quantum state, but only in the 
signs of (cr*), which fully determine the system in the 
classical configuration. This provides some robustness. 
We find that for a successful annealing procedure, yield¬ 
ing the correct sign for all i, the number of phonons pro¬ 
duced during the evolution should not be larger than 1. 
At fixed detuning, we can reduce the number of phonons 
by decreasing the Rabi frequency, at the expense of in¬ 
creasing time scales. As a realistic choice [Til . [TEj ]. we 
demand that annealing is achieved within tens of milisec- 
onds. 

Fig. ©a) shows that one can operate at a de¬ 
tuning S = 2tt x 239 kHz, that is, at the on¬ 
set of a glassy phase according to Fig. [3] The 
mode vector which selects the ground state is £ 5 = 
(0.61,0.34,0.11,—0.11,—0.34,—0.61), and the corre¬ 
sponding pattern can be read out after an annealing time 
t > 15 ms. On the other hand, even for long times, (cr* =4 ) 
saturates only at —0.15, which is far from the classical 
value —1. As shown in Fig. ©b), a slower annealing pro¬ 
tocol leads to more robust results (|(<r* )| > 0.52 V *). In 
Fig. me), a much simpler instance in the ferromagnetic 
regime is considered. Good results (|(cr*)| > 0.36 V i ) 
can then be obtained within only a few ms. 

In addition, dephasing due to instabilities of applied 



t (ms) 


t (ms) 


t (ms) 


FIG. 4: Simulated quantum annealing in the closed 
system. Unitary time evolution of the full system (TV = 
6) for different laser detunings 6 between the 4th and 5th 
resonance. Thick lines show the result of the exact evolution, 
while the thin lines have been obtained from the semi-classical 

approximation. In all cases, the desired pattern (+ + H-) 

can be read out after sufficiently long annealing times. While 
in (a,b), we operate at the onset of glassiness, d = 2-7T x 239 
kHz, the panel (c) considers a ferromagnetic instance, <5 = 
2-7T x 143 kHz. The annealing protocol in (a) is defined by 
B max = 50 J rms and t = 3 ms. In (b), the same instance is 
solved with higher fidelity by choosing R max = 80J rm s and 
t = 6 ms. Fast annealing, with _B max = 50J rm a and r = 1 
ms, is possible in the ferromagnetic instance in (c). In all 
simulations, we have chosen ff = 2-7T x 50 kHz, w reC oii = 27r x 15 
kHz, and e = —10 kHz. 


fields and spontaneous emission processesundesired pro- 
cesses disturb the dynamics of the spins. In Ref. [42| a 
master equation was derived that takes into account such 
noisy environment. To study the evolution of the system 
in this open scenario we have applied the Monte Carlo 
wave-function method [45j . As quantum jump operators 
are hermitean, cr x for dephasing and a\ for spontaneous 
emission, the evolution remains unitary, but is randomly 
interrupted by quantum jumps. Each jump has equal 
probability T, and the average number of jumps within 
the annealing time T is given by Uj umps = 2 NTT, which 
we chose close to 1. 

Since a faithful description requires statistics over 
many runs, we restrict ourselves to a small system, 
N = 4, with short annealing times. In a sample of 100 
runs, we noted 94 jumps (42 cr x and 52 a z jumps). In 39 
runs, no jump occured. Amongst the 61 runs in which at 
least one jump occured, 26 runs still produced the correct 
sign for all spin averages (a x ). The full time evolution, 
averaged over all runs, is shown in Fig. [5] On average, 
the final result is (0.65,0.50,-0.41,-0.68), that is, our 
annealing with noise still produces the correct answer, 
but with lower fidelity. 

Whether an individual jump harms the evolution cru¬ 
cially depends on the time at which it occurs: While a 
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t (ms) 

FIG. 5: Simulated quantum annealing in the open 
system. Using a Monte Carlo wave-function description 
averaged over 100 runs, we simulate the annealing process 
with IV = 4 ions in the presence of dephasing and spon¬ 
taneous emission. With a noise rate F = 0.03J rm s ~ 120 
Hz and a total annealing time T = 1 ms, we have on av¬ 
erage one jump per run, rij umpa = 0.96. We plot expecta¬ 
tion values of the spins and total number of phonons, pro¬ 
duced by the coupling, as a function of time. Here, we have 
operated in the ferromagnetic regime between the second 
and third phonon resonance, S = 2n x 159 kHz. Choosing 
If = 2-7T x 50 kHz, tLi reC oii = 27 t x 15 kHz, u) z = 2n x 876 kHz, 
Umax = 50Jrms ~ 200 kHz, and e = — 10 kHz, we are able to 
perform very fast annealing, t = 0.1 ms, with high fidelity. 
After a time T = 1 ms, the unitary time evolution has con¬ 
verged to ((al),{a 2 x ),(al),{a$)) = (0.97,0.70,-0.65,-0.85), 
correctly reproducing the classical pattern (+ H-). 

spin flip (cr~ noise) is harmless in the beginning of the 
annealing, a dephasing event (a x noise) at an early stage 
of the evolution leads to wrong results. Oppositely, at 
the end of the annealing procedure, dephasing becomes 
harmless while spontaneous emission falsifies the result. 
An optimal annealing protocol has to balance the effect 
of different noise sources against non-adiabatic effects in 
the unitary evolution. 


Scalability 

Above we have demonstrated the feasibility of the pro¬ 
posed quantum annealing scheme in small systems. The 
usefulness of the approach, however, depends crucially 
on its behavior upon increasing the system size. While 
the exact treatment of the dynamics becomes intractable 
for longer chains, an efficient description can be derived 
from the Heisenberg equations: 

ih^( a m) = ([a m ,H(t)]}, (12) 

= {[<7a,H(t)]), 


with H = H 0 (t) + Hs(t) + hbias- To solve this set of 5 N 
first-order differential equations, we make a semi-classical 
approximation (a m <r*) ~ ( a m )(a l x ), and then proceed nu¬ 
merically using a fourth order Runge-Kutta algorithm. 
The semi-classical approximation is justified as long as 
the system remains close to the phonon vacuum. A di¬ 
rect comparison with exact results for six ions (see Fig. 
11 shows that the semiclassical approach, while slightly 
overestimating fidelities, accurately reproduces all rele¬ 
vant time scales. 

This approach allows us to extend our simulations up 
to N = 22 ions, at trap frequency oj z = 2n x 270 kHz. 
We operate between the first and second resonance where 
the level spacing is largest, at a beatnote frequency wl = 
u i + 0.2 (u >2 — uq), that is with a fixed relative detuning 
between the two modes. This choice corresponds to 6 = 
2tt x 1.2 MHz in Fig.d characterized as a glassy instance 
of the system. 

Our aim is to find the relation of annealing time, mea¬ 
sured by the decay parameter r, and systems size while 
the fidelity F is kept constant. For a practical definition 
we demand that F is zero when the annealing fails, that 
is when the sign of the spin averages {a l x ) does not agree 
with the classical target state for all i. If the annealing 
finds the correct signs, the robustness still depends on the 
absolute values of the spin averages. The fidelity is then 
defined as the smallest absolute value, F = ruin, | ) |. 

Our results are summarized in Fig. [G] Firstly, this figure 
shows that for all sizes N < 20 large fidelities F > 0.5 can 
be produced within experimentally feasible time scales, 
r < 30 ms. Secondly, the time scale r needed for a fi¬ 
delity F = 0.5 fits well to a fourth-order polynomial in 
N (with subleading terms of the order exp(l/JV)): 

t(N) = N 4 tq exp( 7 /A r ), (13) 

with To and 7 being free fit parameters. Although the 
sample of 22 ions is too small to draw strong conclu¬ 
sions, it is noteworthy that the polynomial fit is more 
accurate than an exponential one. This suggests that 
the proposed quantum simulation is indeed an efficient 
way of solving a complex computational problem. One 
should also keep in mind that our estimates, based on 
a semi-classical approximation, neglect certain quantum 
fluctuations which could further speed-up the annealing 
process. 

To study the scaling of dissipative effects, we have 
extendend the Monte Carlo wave function approach to 
larger systems, which is feasible if the phonon dynam¬ 
ics is neglected. The unitary part of the evolution is 
then described by the effective Ising Hamiltonian H e g = 
Hj + ffe(t) + hbias- The dissipative part consists of ran¬ 
dom quantum jumps described by <r* and a l z . The results 
for a glassy instance (6 = 2tt x 198 kHz at w z = 2tt x 700 
kHz) are summarized in Table|I]for IV = 4, 6 , 8 . The noise 
rate is chosen such that on average one quantum jump 
occurs in the system with four ions, while accordingly 
the system with eight ions suffers on average from two 
such events. In all cases, the annealing produces the cor- 














FIG. 6 : Fidelity estimates of the quantum annealing 
and scaling of annealing times. We solve the equations 
of motion, Eq. urn for a glassy instance at different system 
sizes TV, and plot the fidelity of the outcome as a function of 
the annealing time t. In the inset, we investigate the scaling 
behavior by plotting (in double-logarithmic scale) the value 
of t which is needed for a fidelity F = 0.5 as a function of TV. 
A fourth-order polynomial fit agrees very well with the data 
(black dashed line). The fit parameters, as defined by Eq. 
(1131 . are to = (90 ± 40) ms and 7 = 12.0 ± 1.2. For all cal¬ 
culations, we have chosen a beatnote frequency between the 
two lowest resonances, col = 0.8u>i + 0 . 2 cj 2 , a trap frequency 
ojz = 2-7T x 270 kHz, and a bias potential e = — 1 kHz. The 
initial value of the transverse field was Umax = 10 kHz. 

rect pattern, F > 0. As expected, F decreases for larger 
systems, but fortunately rather slowly (from F = 0.25 
at TV = 4 to F = 0.16 at TV = 8). If the total amount 
of noise is kept constant, i.e. T oc 1/TV, the annealing is 
found to profit from larger system sizes, since a quantum 
jump at spin i is unlikely to affect the sign of (a J x ) for 
j ^ i. We note that the spin values produced by the 
Monte Carlo wave function method cannot be described 
by a normal distribution. Importantly, the peak of each 
distribution, roughly coinciding with its median, is barely 
affected by the noise. Thus, larger fidelities can be ob¬ 
tained from the median rather than from the arithmetic 
mean of (ex*). 

Spin pattern in the quantum Mattis model 

The quantum annealing discussed above exploits quan¬ 
tum effects in order to extract information encoded in the 
classical model. Now we search for information which is 
encoded in the quantum, but not in the classical model. 
Therefore, we focus on the ferromagnetic Mattis model, 
which in the classical case keeps a binary memory of a 
spin pattern, that is, of TV bits. Our considerations can 



12 3 4 5 6 

/ 


FIG. 7: From binary to real-valued patterns. The spin 
expectation values A i approach the real values of the mode 
vector when a sufficiently strong transverse magnetic field 
B is present. The shown data was obtained from the effective 
spin Hamiltonian for TV = 6 ions in the ferromagnetic regime, 
5 = 2n x 143 kHz. Deviations from the equality A; = £i m are 
smaller than 0.04. 

also be generalized to the Hopfield model p], which mem¬ 
orizes multiple patterns. We will show how quantum ef¬ 
fects can increase the amount of information encoded by 
these models. 

Recall that in the classical case, the spin pattern was 
defined by a resonant mode in terms of the sign of each 
component. In the quantum case, however, one can¬ 
not simply replace classical spins by quantum averages, 
sign((cr*)). Even in a weak transverse field B , this quan¬ 
tity vanishes due to the Zi symmetry, a x —>• — a x . In¬ 
stead, the pattern is reflected by A,; = (Ti| tx* |SI/ 2 ), where 
llki) and ^ 2 ) are the ground and first excited state. For 
small B , we find numerically sign(Ai) = sign(£, m ). For 
large R, the stronger relation A * = t; lrn holds approxi¬ 
mately, see Fig. 0 Thus, the former binary memory has 
become real-valued. 

To show this behavior, we note that for strong B , the 
ground state is fully polarized along z, and the first ex¬ 
cited state is restricted to the TV-dimensional subspace 
with one spin flipped, that is, S z = a l = N — 2. 
Within this subspace the Hamiltonian Hj is given by 
an TV x TV matrix approximately proportional to Jij = 
for i ^ j, and Ju = constant. Here we neglect 
all but the m-th mode close to resonance. 

It is easy to see that the vector is a ground state 
of the matrix —£im£jm, which differs from J\j only by 
the diagonal elements, which approach unity for large TV. 
The first excited state reads |T , 2 ) = J2iLi 6m I*}, where 
|i) denotes the state in which spin i is flipped relatively 
to all others (in the a z basis). This shows that A i « 6m, 
and the small deviations decrease quickly with TV. 

Measuring \ experimentally is possible by full state 
tomography. The absolute value of \ can be obtained 
via a simple a\ measurement. In the limit of strong in¬ 
fields we have \ = [(1 — (tx*))/2] 1 / 2 . 
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closed system: 

(1.0,0.97,-0.39,-0.42) 

N = A 

open system, mean values: 

(0.79,0.70,-0.25,-0.32) 


open system, median values: 

(1.0,0.97,-0.39,-0.42) 


closed system: 

(1.0, 0.99, 0.76, -0.40, -0.64, -0.66) 

N = 6 

open system, mean values. 

(0.82, 0.68, 0.55, -0.23, -0.43, -0.46) 


open system, median values: 

(1.0, 0.98, 0.76, -0.39, -0.64, -0.66) 


closed system: 

(1.0, 0.99, 0.96, 0.58, -0.31, -0.70, -0.75, -0.77) 

N = 8 

open system, mean values: 

(0.83, 0.69, 0.69, 0.40, -0.16, -0.43, -0.55, -0.53) 


open system, median values: 

(1.0, 0.98, 0.94, 0.58, -0.30, -0.69, -0.75, -0.76) 


TABLE I: Annealing in effective spin model with and without noise. We perform quantum annealing (r = 3 ms and 
T = 30 ms) for a glassy instance (<5 = 27t x 198 kHz at uj z — 2n x 700 kHz, e = — 1 kHz) at different system sizes N, using 
an effective spin model description. The results ((cr),},.. ., are shown for the closed-system dynamics and for a noisy 

system, with mean and median over a sample of 2000 runs. The average number of noisy event scales with the system size, 
and is adjusted to A/N. In all cases, the signs of (crj,) reproduce correctly the mode pattern, and the fidelity decreases with the 
system size. In contrast to the arithmetic mean values, the median values in the noisy sample are barely affected by the noise. 


Many applications are known for the classical spin sys¬ 
tem with couplings defined by spin patterns, reaching 
from pattern recognition and associative memory in the 
Hopfield model [l| to noise-free coding @,[46j]. Our analy¬ 
sis suggests that patterns given by real numbers could re¬ 
place patterns of binary variables by exploiting the quan¬ 
tum character of the spins. 

Discussion 

In summary, our work demonstrates the occurrence of 
Mattis glass behavior in spin chains of trapped ions, if 
the detuning of the spin-phonon coupling is chosen be¬ 
tween two resonances. In these regimes, the effective spin 
system has an exponentially large number of low-energy 
states, and finding its ground state corresponds to solv¬ 
ing a number-partitioning problem. This establishes a 
direct connection between the properties of a physical 
system and the solution of a potentially NP-hard prob¬ 
lem of computer science. Given the state-of-art in ex¬ 
periments with trapped ions, the physical implementa¬ 
tion is feasible: In comparison to previous experiments 
with trapped ions llfTll4j , only the detuning of the spin- 
phonon coupling needs to be adjusted. Differently from 
other approaches to spin glass physics, our scheme does 
not require any disorder. In its most natural implemen¬ 
tation, parity symmetry allows one to analytically de¬ 
termine the ground state. Different ways to break this 
symmetry can be implemented to increase the complexity 
of the problem. 

The ion chain then becomes an ideal test ground for ap¬ 
plying quantum simulation strategies to solve computa¬ 
tionally complex problems. By applying a transverse field 
to the ions, quantum annealing from a paramagnet to the 
glassy ground state is possible. The ionic system may be 
used to benchmark quantum annealing, which has be¬ 
come a subject of very lively and controversial debate 
since the launch of the D-Wave computers j3, d] • Exact 


calculations for small systems (iV = 6) and approxima¬ 
tive calculations for larger system (N = 22) demonstrate 
the feasibility of the proposed quantum annealing, and 
suggest a polynomial scaling of the annealing time. Ac¬ 
cordingly, this approach may offer the sought-after quan¬ 
tum speed-up. In view of sizes of 30 and more ions al¬ 
ready trapped in recent experiments (cf. Ref. [13), a 
realization of our proposal could not only confirm our 
semi-classical results, but also go beyond the sizes con¬ 
sidered here. 

Finally, resonant coupling to multiple modes opens an 
avenue to neural network models, where a finite number 
of patterns is memorized by couplings to different phonon 
modes. Quantum features can increase the memory of 
such networks from binary to real-valued numbers. It 
will be subject of future studies to work out the possible 
benefits which quantum neural networks may establish 
for information processing purposes. 
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